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ABSTRACT 

In this article we investigate the outer and inner mass distributions of the irregular 
galaxies UGC 4284 and UGC 11861, taking advantage of published HI and Ha high 
resolution rotation curves and constraining the stellar disk of both galaxies throughout 
stellar population synthesis studies. In addition we take into account the gas content 
of both galaxies deriving the HI+He rotation curve. The deduced baryonic rotation 
curves (star+gas) are inadequate to account for the total mass of UGC 4284 and 
UGC 11861, for that reason we examine the possibility of dark matter to explain the 
incongruity between the observed HI and Ha rotation curves of UGC 4284 and UGC 
11861 and the derived baryonic rotation curves. We consider NFW, Burkert, DiCintio, 
Einasto, and the Stadel dark matter halos, to analyse the dark matter content of UGC 
4284 and UGC 11861. The principal results of this work are that cored dark matter 
models better reproduce the dark matter Ha and HI rotation curves of UGC 11861 and 
the dark matter HI rotation curve of UGC 4284, while, the Ha rotation curve of UGC 
4284 is better reproduced by a cuspy DiCintio DM model. In general, cored exponential 
two-parameters models Einasto and Stadel, give better fits than Burkert. This trend, 
as well as to confirm past results, presents for the first time a comparison between 
two different exponential dark matter models, Einasto and Stadel, in an attempt to 
better constrain the range of possible exponential dark matter models applied to real 
galaxies. 

Key words: galaxies:irregular — galaxies:kinematics and dynamics - galax- 
ies:individual:UGC 4284 - galaxies:individual:UGC 11861 - (cosmology:) dark matter 


1 INTRODUCTION 


The rotation of gas and/or stars in galaxies provides an 
important information about the kinematics of the galactic 
disk. Observed RCs are somewhat flat and no Keplerian 
decline was detected. The interpretation based on the dark 
matter (DM) paradigm assumes a large amount of DM far 
beyond the visible part of galaxies. |Freeman| ( |1970[ ) was one 
of the first to suggest that RCs could be used as indication 
of DM, reporting that the RCs of NGC 300 and M33 did 
not show the expected Keplerian decline, and estimating 
that the DM mass in these galaxies should be at least as 


large as the mass of the luminous matter. |Rubin et al.| 
(19801 and Bosma ( 1978| 1981a|b ( endorse the evidence of 
DM, performing extensive studies of galaxies RCs in optical 
and radio wavelengths. [Bahcall et al.| ( |1982[ ) built artificial 
RCs considering a central component, a disk, a spheroidal 
component and a DM component to determine unexpected 
deviations from the generally believed logarithmic slope 
of 2.0. Carignan & Freeman (19851 analysed the mass 


distribution of four galaxies detecting the necessity of a DM 
component to explain the observed RCs of those galaxies 
whose RCs extend to 1.5 — 3.0i?25. These works turn out 
a general acceptance of the DM hypothesis. 
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In general, the first step to determine the DM distribu¬ 
tion in galaxies is to recognize what kind (bulge, disk, and 
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bar) and how many stellar components constitute the inner¬ 
most optical disk. This can be done by performing isopho- 
tal analysis on broad band image photometry of the stud¬ 
ied objects. In fact, the luminosity profile obtained from 
galaxies observations, is the nnique information to derive a 
coarse estimate of the visible mass of such objects, assum¬ 
ing a constant relation between luminosity (L) and mass 
(M). Once all the stellar components are identified, differ¬ 
ent mass profiles for the stellar disk, bulge, and bar are 
used to account for the stellar mass of the optical disk. 
The gas mass content is represented by the mass profile 
of the gaseous tracer (usually neutral hydrogen) used to 
measure the RC far from the optical disk extent. The DM 
halo is commonly considered spherically symmetric and its 
circular velocity is given by the ordinary Newtonian law. 
The next step is to accomplish a non-linear least squares 
ht to the observed RC considering a general fitting rela¬ 
tion VtltiR) = [V.RR) + V,IRR) + V^,i,{R)]. The prob¬ 
lem arise for the luminous components, because in general 
the stellar mass is unknown and a guess of the mass-to- 
luminosity ratio (M/L) is required. In contrast, the gas 


amount is usually well-known and rather small (Swaters 
et al.|2000 1. During the fitting procedure, generally the M/L 


is a free and constant parameter given the ignorance about 
the stellar content of a disk galaxy. 

The maximum disk hypothesis (MDH) was introduced 
by van Albada et al. (19851 and van Albada & Sancisi (19861 


to study the DM distribution of NGC 3198 and it is based on 
the assumption that the DM halo dominates at radii greater 
than 2.2 Rn- Inside this radius the luminous matter could 
account for the galaxy circular velocity estimated up to that 


radius. Ratnam & Salucci (20001 analysed the mass distri¬ 


bution of 83 spiral galaxies, using high resolution RCs data, 
formulated a more modern version of the MDH, suggesting 
a baryonic scale range in radius from 0.5 Rd ~ 0.5 kpc to 
2.0 Rd ~ 30.0 kpc, where the luminous matter seems to 
produce a substantial part of the gravitational potential of 
the studied galaxies. 

In general there are two main factors that govern the 
degeneracy of the results of the RCs fitting: actual errors 
in RCs measurements and our ignorance of the true distri¬ 
bution of luminous matter. Hypothetically, htting different 
mass components to observed RCs should be a valid ap¬ 
proach if we have a great number of high resolution data 
with negligible measurement errors. In practice the observed 
RCs are far from these ideal conditions, therefore the RCs 
decomposition continue to be a degenerate problem. A more 
robust determination of the distribution of the luminous 
matter based on more accurate constraints of the M/L vari¬ 
ation with galactocentric radius and stellar colours could 
lessen the degree of degeneracy. 

Some authors found a strong correlation between the 
variation of stellar M/L in optical and near-infrared (NIR) 


passbands and stellar population colours (Bell & de Jong 
2001 Bell et al.|2003 1. Bell & de Jong (20011 also discovered 


that low luminosity galaxies have lower stellar M/L than 
high surface brightness galaxies and galaxies with few gas 
content. From these works it remains clear that to remove 
the above mentioned degeneracy, due to our poor knowl¬ 
edge of the stellar M/L, more information about the radial 
variation of the stellar M/L and the correlation with other 
galaxies properties are necessary. Gallazzi et al. (20051 ex¬ 


plore the dependency among the total galaxy stellar mass, 
the stellar metallicity, and the stellar age, estimating these 
parameters for a magnitude-limited sample of ~ 200,000 
galaxies selected from SDSS DR2. These authors made a 
comparison of the observed galaxy spectra and a large li¬ 
brary of synthetic spectra based on Monte Carlo library of 


star formation histories (SFHs). Gallazzi et al. (2005 I and 


Gallazzi et al. (20081 found a large scatter in the relations 


amidst stellar mass, age and metallicity, concluding that the 
physical parameters of galaxies at the present time are not 
uniquely determined by stellar masses, and that gas infall or 
outflow could also play an important role in building galactic 
physical properties. 

From these works it is evident that the stellar M/L rely 
on so many physical parameters that the simplistic assump¬ 
tion of considering the stellar M/L constant along the galac¬ 
tic disk is unreal. In this respect, stellar population synthesis 
models (SPs) represent a valuable tool to constrain in a more 
realistic manner the true value of the stellar M/L. |Zibetti| 


et al. (20091 conceive a novel approach to build spatially 


resolved maps of stellar mass surface density considering 
broad-band photometry images of galaxies in the optical and 
NIR. The authors used a Monte Garlo library of 50,000 syn¬ 


thetic SPs by Bruzual & Chariot (20031 and Bruzual (20071 


incorporating an updated recipe for the thermally pulsing 
asymptotic giant branch stellar evolutionary phase following 


Marigo & Girardi (20071 and Marigo et al. (20081. Gallazzi 
fc Belli (|2009[) studied the accuracy in the determination 


of stellar M/L from galaxy colours, using the SPs models 


of Bruzual & Chariot (2003 \ and a Monte Carlo library of 


150,000 SFHs, to examine which spectroscopic conditions 
and spectral signal-to-noise ratios (S/N) are required to de¬ 
crease the errors on M/L determination for different galaxy 
types. The authors concluded that in galaxies with old stel¬ 
lar populations the accuracy of M/L is beneath 0.05 dex for 
S/N > 20, while for galaxies with young stellar population 
the M/L presents, in general, higher statistical uncertain¬ 
ties. The results of Gallazzi & Bell (20091 are based on the 


very strong assumption that the dust extinction is negligible. 
|MacArthur| ( |2005[ ) prove that the effect of dust attenuation 
is important near 4000 A, additionally, from the literature 
it is well-known that dust effects have a crucial impact on 
M/L values based on colours ( [Zibetti et al.|[2C)09| . For this 
reason, at the moment, it does not exist a robust analy¬ 
sis discussing the feasibility of SPs studies to determine the 
M/L, nonetheless SPs seem to be the more powerful tool to 
obtain the stellar mass content of galaxies. 


The A-CDM paradigm (Blumenthal et al. 19841 predicted 
cuspy DM density profiles in the centre of galaxies ( [Navarro] 
et al.|1996 1 and recent results with higher resolution seem to 


suggest a shallower DM density profile (Navarro et al.|2004 


|2010p . Nonetheless the past and newer developments are in 
contradiction with observational results of cored DM den¬ 
sity profiles in dwarf galaxies (de Blok et al.||2001 de Blok 


&Bosma2002 Swaters et al. 2003 Ghemin et al. 20111. This 


discrepancy is known as the cuspy-core problem. Many other 
authors, in an attempt to contribute to the solution of the 
cuspy-core problem, performed the RG decomposition con- 


sidering a constant M/L ( 

Carignan & Freeman 

1985 

1988 

Jobin & Carignan|19901 [Persic & Salucci|1990 [Martimbeau 

et af.|1994 

Persic et af.|1996| Bfais-Oueffette et af.|2001| de 

Blok et al. 

2001| de Blok & Bosma|2002| Swaters et al.|2003| 
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Kuzio 

de Naray et 

al.||2006 [Fuentes-Carrera et al.||2007| 

Spano et al.|2008| Kuzio de Naray et al.|2008||Repetto et al.j 

2010 

Chemin et al. 

2011 

1 and reinforce the evidence of a 


cored DM distribution in the central part of galaxies. Con¬ 
versely, other authors pursue spectrophotometric and SPs 


studies to derive the M/L of the stellar component (Salucci 
et al. 12008) [Repetto et al.|2013| |de Denus-Baillargeon et al. 


2013) or at least to constrain the stellar component (Dut- 
ton et ar]|2005 I to avoid the disk-halo degeneracy. Salucci 


et al. (20081 accomplished RCs mass modelling and spectral 


energy distribution fitting with spectrophotometric models 
to obtain the disk masses of 18 Sa spiral galaxies (princi¬ 
pally bulgeless systems) finding that by decomposing the 
RCs with the spectrophotometric disk masses the results 
are consistent with the corresponding maximum disk solu¬ 


tion. Repetto et al. (20131 employed the Zibetti et al. (20091 


method to optical (SDSS) and NIR (2MASS) images of NGC 
5278 (KPG 390A) to obtain the stellar disk mass profile of 
that galaxy from SPs models to reproduce the RC of KPG 
390A. The new strategy relied on fitting only the DM RC, 
obtained by subtracting the SPs baryonic disk from the ob¬ 
served RC of NGC 5278. The most important finding of 


Repetto et al. (2013) is that the favoured DM distribution 


is cored when the disk mass approximate to the maximum 
disk solution in agreement with the general belief, [de Penns-] 


Baillargeon et al. (2013) used a chemo-spectrophotometric 


galactic evolution model to determine the stellar M/L and 
perform the RC decomposition of ten spiral and dwarf ir¬ 


regular galaxies from the SINGS survey ( Kennicutt et al. 
|2003[ ). The authors employed the settled baryonic disk as a 
weighting function to fit the model DM halo to the RCs of 
the studied subsample of galaxies, and concluded that the 
stellar disks obtained from their chemo-spectrophotometric 
models were compatible with the maximum disk hypothesis. 
The few examples presented, indicate that a growing effort 
to break the disk-halo degeneracy exists, however the most 
signihcant contributions trying to address the cuspy-core 
problem still rely on the general assumption of considering 
the stellar M/L constant along the galactic disk. In general, 
it is still missing a significant endeavour to earn the DM dis¬ 
tribution, through RC fitting, determining the disk stellar 
mass from SPs studies. For this reason it is worth consid¬ 
ering the cuspy-core problem with a different observational 
approach, focusing on the formulation of a general proce¬ 
dure to better constrain the baryonic disk mass in galaxies. 
This is the first of a series of papers dedicated to the study 
of the outer and inner mass distribution in irregular galax¬ 
ies, applying and testing the novel method of |Repetto et ah] 
120131, using the RCs data of the WHISP HI survey of galax¬ 
ies ( Verheijen fc Sancisi]2001 1, to investigate the outer mass 
distribution, and the GHASP Ha survey of galaxies, to anal- 


yse the inner mass distribution 

Amram et al. 

2002 

Garrido 

et al. 2002 Spano et al. 2008 

Epinat et al. 

2008 

1. In this 


first article we concentrate on the determination of the mass 
distribution of two irregular galaxies UGC 4284 and UGC 
11861, whose RCs are somewhat symmetric and regular and 
for this reason more suitable to perform the RC analysis. 
The novelty of the approach conceived in [Repetto et al.| 
(2013) resides in deriving the stellar disk mass of the galaxy 
from SPs studies and in subtracting the built baryonic disk 
from the observed RC to fit only the remaining DM RC as 
we have already summarised above. 


The main motivation of these works are attempt to answer 
the following question: Performing the RC decomposition 
with more realistic M/L obtained from SPs studies, could 
change the determination of the DM halo parameters sub¬ 
stantially? In the case of a positive answer, this research 
could shed new light on how the DM halo parameters change 
because of a non parametric baryonic component. In the case 
of a negative answer the analysis is still worth to characterise 
a novel procedure to constrain the stellar disk parameters 
with a deeper knowledge of the physics of the underlying 
stellar disk populations derived from photometric data. 

This paper is organized as follows. Sections 2 and 3 
describe the GALFIT disk/bulge/bar decomposition. Sec¬ 
tion 4 describes the constrains of the baryonic M/L ratio 
through SPs. Section 5 describes the Ha and HI RC deriva¬ 
tions. Section 6 presents the analysis of the inner and outer 
DM distributions. Section 7 discusses the results. Finally, 
Section 8 presents our conclusions. 


2 2D GALFIT DISK/BULGE/BAR 
DECOMPOSITION OF UGC 4284 


UGC 4284 is classified as Scd(s) in the RC3 (de Vau- 
couleurs et al.||1991|), and as SA Bc by|Paturel et al.|(|2003p . 
Kandalian & Kalloghlian ( 1998|), Chapelon et al. (19991 and 


Cabrera-Lavers & Garzon (20041 found evidence of a cen¬ 


tral isophotal distortion that they explain postulating the 
existence of a bar with diameter of 38”. Other authors did 
not confirm the presence of a bar in UGC 4284, nevertheless 
they encountered non-circular and off-plane gas motion in 
the inner region of UGC 4284 that they discuss as originated 
by a warped gaseous disk ( Coccato et al.|2()04 J6zsa|20W I. 
Bearing in mind the purpose of the present study it is im¬ 
portant to determine which kind of baryonic components do 
actually exist for UGC 4284, and given the incongruity of the 
above reported results we required to perform 2D isophotal 
analysis of the surface brightness of this galaxy, using GAL¬ 
FIT ( Peng et al.|2002 2010| and considering at least three 
principal components: disk, bulge and bar. UGC 4284 has 
been observed with Hubble Space Telescope (HST) in the 
filters F555W (WFPC2 V), F814W (WFPC2 wide I) and 
F606W (WFPC2 wide V), nonetheless the filters F555W 
and F814W show only a small part of the periphery of the 
galaxy. In contrast the filter F606W exhibits a larger portion 
of the disk of UGC 4284, so we use that filter to decompose 
the disk of UGC 4284. The results of the 2D exponential disk 
fit are presented in Table 1. In the case of UGC 4284 we also 
examined the possibility of a bulge and a bar components 
given the inconsistency in the morphological classification of 
this galaxy. We use the Spitzer IRAC 4.5 pm image because 
UGC 4284 is indiscernible in the 2MASS J, H, K bands, and 
in all the SDSS filters is too faint to allow an obvious sep¬ 
aration from the background. Neither the Sersic component 


nor the Ferrers | Ferrers 18771 component produce accept¬ 


able fit parameters to reveal the presence of a bulge and a 
bar, for this reason we conclude that UGC 4284 has only 
an exponential disk component in total agreement with the 
RC3 classification. 
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4 P. Repetto, et al. 


Table 1. Results of GALFIT isophotal analysis for UGC 4284. 


Table 2. Results of GALFIT isophotal analysis for UGG 11861. 


UGC 4284 exponential disk final parameters. 

Surface Brightness (mag arcsec“^) 

17.97 

Scale length (kpc) 

1.27 

Axial ratio 

0.46 

Inclination (degrees) 

64.98 

Position Angle (degrees) 

44.45 


UGC 11861 exponential disk final parameters. 

Surface Brightness (mag arcsec“^) 

19.35 

Scale length (kpc) 

1.30 

Axial ratio 

0.47 

Inclination (degrees) 

64.27 

Position Angle (degrees) 

71.46 


UGC 11861 Ferrers bar final parameters. 


3 2D GALFIT DISK/BULGE/BAR 
DEGOMPOSITION OF UGG 11861 


UGC 11861 was classified as SABdm in the RC3 |de 
|Vaucouleurs et al.|.|1991| , as SABd by |Paturel et al.] ( |2003[ ). 
At first glance the presence of a bar is not so evident from 
single optical (SDSS) and NIR (2MASS) broad band im¬ 
ages, nevertheless a central boxy structnre is clearly visible 
in the Spitzer IRAC 3.6 ^m image. The congruity in the 
morphological taxonomy of UGC 11861 allow us to identify 
promptly at least two components for that galaxy: a disk 
and a bar. We use GALFIT to prove the existence of such 
structures for UGC 11861. We employ the SDSS g band im¬ 
age for the disk and the Spitzer IRAC 3.6 fj,m image for 
the bar. We consider an exponential disk and Ferrers model 
( Ferrers|[l877 1 for the bar. The GALFIT result for the disk 
and bar GALFIT model are listed in Table 2. We also con¬ 
sider the possibility of a bulge component for UGC 11861, 
for which we adopt a Sersic profile and we fit with GALFIT 
2MASS J, H, and K band images regarding a de Vaucouleurs 
bulge with an older stellar population, and SDSS g,r and i 
band images considering the possible presence of pseudob¬ 
ulges with younger stellar population. The choice of that 
photometric range to examine the possible existence of a 
bulge for UGC 11861 is also dictated by the fact that the 
bar emission is fainter in the selected passbands, and this 
fact could in principle allow an unambiguous detection of a 
potential bulge. In the case of the SDSS g, r and i bands we 
first fit an exponential disk model and then we try to fit a 
Sersic component to detect the existence of an hypothetical 
bulge in the residuals without success. The GALFIT fit to 
the 2MASS J, H, and K passbands did not reveal the pres¬ 
ence of a bulge. From the GALFIT analysis we conclude 
that UGC 11861 is composed of an exponential disk and a 
bar, in good agreement with the morphological classification 
of that galaxy. 


4 BARYONIC M/L RATIO FROM RESOLVED 
MASS MAPS 


Reliable M/L ratios are an important aspect when de¬ 
termining the baryonic contribution to the rotation curves of 
disk galaxies. By using a chemo-spectrophotometric galac¬ 
tic evolution model [de Denus-Baillargeon et al.| | |2013P high¬ 
lighted the importance of adopting a radially varying M/L 


ratio, vs. a radially constant one (see also Repetto et al. 


In this work we adopt a resolved mass map method (Zi- 


betti et al. 20091 that relies on a Monte Carlo library 


Surface Brightness (mag arcsec“^) 

8.83 

Outer truncation radius (kpc) 

2.48 

a. 

2.95 

/3 

0.0 

Axial ratio 

0.38 

Inclination (degrees) 

70.74 

Position Angle (degrees) 

-16.6 


of 50,000 stellar population spectra randomly drawn from 
the SPS models of S. Chariot & G. Bruzual (CB07). 
These models use single stellar populations as described 
by Bruzual & Chariot 120031, but include the prescription 
for stellar evolution in the thermally-pulsating asymptotic 
giant branch (TP-AGB) by Marigo & Girardi (20071. The 
IMF of Chabrier (Galaxy-disk, 20031 is adopted. The ad¬ 
justable parameters of the models are the star formation 
history (SFH), which is exponentially declining with ran¬ 
dom bursts superimposed, the metallicity, uniformly dis¬ 
tributed between 0.02 and 2 times solar, and the extinction 
by dust, b ased on the two-component model of [Chariot fc| 


Fall 


(2000| ^The Monte Carlo library is obtained for a given 


combination of colours, and the mass-to-light ratio (M/L) 
for a given filter. 


Zibetti et al. (20091 fiducial method involves the use of 


the g and i bands from SDSS, and one NIR filter. The main 
reason for the use of this combination of filters is that it pro¬ 
vides the less degeneration for the SPs models in a colour- 
colour diagram (e.g., g — i, ys. i — H). Other combinations 
of filters, e.g, Johnson’s V and I combined with one NIR fil¬ 
ter, show a smaller area in the colour-colour diagrams, thus 
providing more degeneracy, i.e., the same colours can be at¬ 
tributed to many different model parameters. The same ap¬ 
plies to various other combinations of filters (see figure]^. It 
was found that the absolute mass values are overestimated, 
as compared to the fiducial method, when a combination of 
colours with more degeneracy is adopted. 

The M/L is obtained for the NIR filter as a function 
of (g — i) and, e.g., (i — H). The models are grouped in 
bins that take into account the observational error in the 
photometry, typically 0.05 mag. For each bin, the median 
value of M/L is computed. The final look-up tables can be 
compared with observations in a pixel-by-pixel basis. A map 
of M/L is obtained which can be converted to absolute mass 
if the distance to the object is known. 

When analysing the spectral energy distribution (SED) 
of galaxies, the wavelength range between 2.5p.m and lOpm 


^ The model takes into account the fact that stars are born 
in dense molecular clouds. A different treatment is given for the 
attenuations of young and old stars. 
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is considered to be a mixture of stellar and dust emis¬ 


sion (e.g., da Cunha et al. 20081. In order to apply the 
resolved mass map method (Zibetti et al. 20091 to get a 


M/L map of, e.g., the IRAC 3.6p.m filter, an infrared SED 


model is needed. The model of da Cunha et al. (20081 is ad¬ 


equate for this purpose. This model incorporates four main 
dust components, the emission from polycyclic aromatic hy¬ 
drocarbons (PAHs), the mid-IR continuum from hot dust, 
and both the emission from warm and cold dust in thermal 
equilibrium. An energy balance condition is applied in a way 
that the total energy absorbed by dust, in stellar birth clouds 
and in the diffuse ISM, is re-radiated in wavelengths over the 
range from 3p.m to 1000p.m. It is assumed that the dust in 
the ISM is heated only by starlight. The infrared model can 
be easily combined with the SPs models of |Bruzual fc Char-| 
lot (20031 and CB07. A Monte Carlo library of 50,000 dust 


emission spectra, for wide ranges of dust temperatures and 
fractional contributions by the different d ust components, 
is provided in the MAGPHY^ package of da Cunha et al. 


(20081. We combine this library with the SPs library tak¬ 


ing care of the energy balance condition as described in |da| 


Cunha et al.|(2008|. A total of 667 million spectra were ob¬ 


tained. The models were grouped in bins of 0.05 x 0.05 mag , 
for the (g — i), and (i —3.6pm) colours. The mean M/L ratio 
for the 3.6pm filter was obtained for each bin. The results 
are shown in figure]^ together with the rms values. 

We then use these results (in figure to obtain the 
two-dimensional stellar mass map, for UGC 11861. In hg- 
ure 1^ we show the azimuthally-averaged radial profiles of 
the M/L( 3 , 6 pm) ratio (left plot), and the stellar (p,) mass 
surface density (right plot). We assume an object distance 
of 18.4 Mpc, and correct the models for Galactic extinc¬ 


tion ( Schlafly fc Finkbeiner|2011| Chapman et al.||2009 K 
For the case of UGC 4284 the resolved mass map 
method was also adopted. The 3.6pm data was taken from 


the Spitzer Survey of Stellar Structures in galaxies (Sheth 
|et al.|[^10[ ). The resulting azimuthally averaged radial pro- 
Hies for the M/La.epm ratio, and the stellar mass surface 
density (p*) are shown in figure]^ We assume a distance of 
9.42 Mpc. 


4.1 Baryonic HI-fHe gas components of UGC 
4284 and UGC 11861 


The HI gaseous RCs of both galaxies were obtained in¬ 
tegrating the total HI column density map in elliptical rings 
using the task ELLINT in the GIPSY package. The outputs 
are the HI mass surface density and the HI cumulative mass 
displayed in figures and The contribution of He and 
metals was taken into account scaling by a factor of 1.4 the 
HI surface density and mass. The total HI-|-He RCs of both 
galaxies are displayed in figure 


5 Ha AND HI RC OF UGC 4284 AND UGC 
11861 


^ http://www.iap.fr/magphys/magphys/MAGPHYS.html 
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Figure 1. Area covered in a colour-colour diagram, for different 
colour combinations obtained from a Monte Carlo SPs library 
grouped in colour bins of 0.05 X 0.05 mag^, as described in |ZibettI| 
|et al.| (|2009|l. Filters u, g, i, and z from SDSS. Filters J, H, and 
K from 2MASS. 


The Ha RCs of UGC 4284 and UGC 11861 were derived 
in the frame of the Gassendi Ha survey of spirals (GHASP 


Garrido et al. 20021, dedicated to study the 2D Ha veloc¬ 


ity field and RCs of a sample of 203 spiral and irregular 
galaxies. [Garrido et al.| ( |2003| |2004| |2005| ) perform the first 
data reduction and analysis of 97 galaxies in the GHASP 
sample. Epinat et al. (20081 accomplished again the reduc¬ 
tion process of the 97 galaxies more 108 new galaxies, using 
adaptive binning technique based on Voronoi tessellation in 
order to provide different weights for low and high signal-to- 
noise zones in the velocity fields. The final high resolution 
2D velocity fields of the GHASP survey have a spectral ac¬ 
curacy of about 5 km s“^ and a spatial resolution of ~ 2”. 
From those high resolution 2D velocity fields, [Epinat et al.j 
(20081 constructed the RCs for the most galaxies of that 
survey. A part of the GHASP sample (130 galaxies) has 
been previously observed by the Westerbork survey of HI in 
spiral galaxies (WHISP) intended to investigate the 2D HI 


kinematic of 500 spiral and irregular galaxies (Verheijen & 
[Sancisi|2001[ ). The typical spectral resolution of the WHISP 
sample of 2D HI velocity fields is about 5 km s“^ with a 
spatial resolution of ~ 12". 

[van Eymeren et al.[ ( [2011[ ) obtained the HI RGs of a sam¬ 
ple of 70 spiral and irregular galaxies, also including UGC 
4284 and UGC 11861, using the WHISP HI velocity fields 
by applying a tilted-ring analysis to derive accurately the 
kinematic parameters from the velocity fields, obtaining 
residuals of a few km s“^. In the particular case of UGC 
4284 and UGC 11861, the HI WHISP RCs are relatively 
more extended that the Ha RCs, they are regular and in 
the inner regions they have a good accordance with the 
Ha RC. We decide to use the Ha RCs to explore the in¬ 
ner DM distribution and the HI RCs to study the outer 
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Figure 2. Left. Logarithm of the effective, i.e., as seen by the observer, M/L ratio for the 3.6|rm filter, M/L (3 in a (g — i), vs., 

(i — 3.6(rm), colour-colour diagram. Data are grouped in bins of 0.05 X 0.05 mag^. Right. Logarithm of the root mean square (rms) errors, 
corresponding to the M/L (3 data in the left panel. The mean error value (of all plotted data) is ~ 30%. 




Figure 3. Radial profiles for UGC 4284. Left Azimuthal average of the stellar mass surface density, as a function of radius. Right 
Azimuthal average of the mass-to-light ratio, M/Ag gi^mi a function of radius. An inclination angle of 59°, and position angle of 
169.9°, were adopted as deprojection parameters. 


DM distribution. The comparison between the Ha and HI 
RCs is shown in figure In the same figure are also pre¬ 
sented the baryonic RCs, considering errors giving a varia¬ 
tion of more and less 30% of the disk mass for both galax¬ 
ies. From figure it is clear that stars and gas alone are 
not enough to account for the total mass of UGC 4284 


and UGC 11861. The total disk mass of UGC 4284 is 
Md = 7.92 X 10®Mq, the total disk mass of UGC 4284 more 
30% is Md = 1.03 X IO^^Mq and the total disk mass of UGC 
4284 less 30% is Md = 5.55 x 10 ®Mq. The total disk mass of 
UGC 11861 is Md = 5.38 x 1O^®M0, the total disk mass of 
UGC 11861 more 30% is Md = 6.99 x 1O^®M0 and the total 
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Figure 4. Radial profiles for UGC 11861. Left: Azimuthal average of the stellar mass surface density, /r*, as a function of radius. Right. 
Azimuthal average of the mass-to-light ratio, M/L (3 as a function of radius. An inclination angle of 61°, and position angle of 

10.7°, were adopted as deprojection parameters. 



Figure 5. HI mass gas profiles for UGC 4284. Left. Azimuthal average of the HI mass surface density, T,uj, as a function of radius. 
Right. Azimuthal average of the HI mass, Mjji (10^ Af©), as a function of radius. An inclination angle of 59°, and position angle of 
169.9°, were adopted as deprojection parameters. 


disk mass of UGC 11861 less 30% is Md = 3.76 x lO^^M©. 
In the case of UGC 11861 the disk mass is the sum of the 
disk and the bar component (see section]^. In the next two 
sections we explore the possibility of DM, for UGC 4284 
and UGC 11861, to explain the discrepancy between the to¬ 


tal observed mass distribution and the mass content due to 
stars and gas. 
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Figure 6. HI mass gas profiles for UGC 11861. Left: Azimuthal average of the HI mass surface density, as a function of radius. 

Right: Azimuthal average of the HI mass, Mjfi (10® Mq), as a function of radius. An inclination angle of 61°, and position angle of 
10.7°, were adopted as deprojection parameters. 


6 INNER AND OUTER DM DISTRIBUTION 
FOR UGC 4284 AND UGC 11861 


We investigate the inner and outer DM content of UGC 
4284 and UGC 11861 accomplishing the Ha and the HI DM 
RCs fitting for both galaxies, considering possible errors in 
the stellar disk determination of plus and minus 30% for 
all the DM models examined in this study. The HI DM 
RCs were derived by subtracting the stellar RCs and also 
the HI+He gaseous RCs from the observed HI RCs of both 
galaxies. 

We consider the DM halo of [Navarro et ^ ( |1996[ ) 
(NEW), the Burkert (BK) Burkert ( 1995[ ) DM halo and the 


Di Cintio et al. (20141 DM halo (DC). The DC DM halo 


was formulated from the double power-law model of |Zhao| 
119961 expressing the three power-law exponents as func¬ 
tional forms of the stellar-to-halo mass ratio, establishing a 
two free parameters DM halo. We use the prescription of 


Dutton & Maccio (20141 for the concentration parameter 


for the NEW DM halo. 

In the present study we use the minuit minimizer 
through the ROOT package written in C++. The ROOT 
package is a huge object-oriented data analysis frame¬ 
work developed at CERN that allows sophisticated anal¬ 
ysis of practically any kind of data. The minuit mini¬ 
mizer uses the MIGRAD algorithm that implements a vari¬ 
able metric method, with inexact line search checking for 
positive-definiteness. The main difference with respect to 
the Levemberg-Marquardt method is that, the minimization 
algorithm only converges on the global minimum, for this 
purpose an important output parameter is the estimated 
distance from the (global) minimum (Edm), a quantity that 
informs the user when the global minimum is reached. The 


Edm together with the values and also the filled con¬ 
tours, represent the fundamental criteria used in this study 
to assess the reliance of the determined solutions. 

The fitting procedure employed in this study considers 
as free parameters the DM halos masses and radii respec¬ 
tively within the intervals [10®, 10^^] Mq and [0.1, 100.0] 
kpc. The solutions we have found are listed in tables 2, 
3, 4 and 5. In addition, the results with stellar disk mass 
Md = 7.92 X 10® Mq for UGC 4284 and stellar disk mass 
Md = 5.38 X 10^'^ Mq for UGC 11861 are displayed in fig¬ 
ures [^13[ The X® filled contours of the solutions presented 
in figures to [^are exhibited in figures and |15| 

Erom the tables and figures presented in this section it 
is evident that the NEW DM halo does not fit either the 
Ha or the HI DM RCs of UGC 4284 and UGC 11861. The 
Edm indicates that the solution is far from the expected 
global minimum and this fact is also confirmed by the non 
existence of the x® filled contours. The BK DM halo fits the 
Ha and HI DM RCs of UGC 4284 and UGC 11861 better 
than NEW, even if the x^ filled contours are narrow and 
for the HI DM RC of UGC 4284 the final halo mass value 
coincide with the lower halo mass limit allowed by our fitting 
procedure. 

The DC DM halo seems to perform better than BK DM 
halo as indicated by the x^ values, the x^ filled contours and 
the Edm values, even if in the case of the DM RC of UGC 
11861 the x^ values for the BK DM halo are lower than the 
DC DM halo. The DM masses originated from the fitting 
of the DC DM halo are in general a factor 10® — 10® times 
greater than the DM halo masses generated by the fitting 
of the BK DM halo. In the case of the HI DM RC of UGC 
11861 the DM masses of the DC DM halo are a factor of 10'* 
greater than BK DM halo, even if in this case the BK DM 
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Figure 7. Ho and HI RCs of UGC 4284 and UGC 11861 together with the corresponding baryonic stellar and gaseous (HI+He) RCs. 
The comparison between the two tracers is to motivate our choice to use the Ha RCs to analyse the inner DM distribution and the HI 
RCs to estimate the outer DM distribution. In this figure is also shown the preliminary fit to the Ha and HI RCs with the corresponding 
values and residuals. The functional forms used for the preliminary fit are exponential-like for the Ha RCs of UGC 4284 and UGC 
11861 and polynomials for the HI RCs of both galaxies. We used those preliminary fits instead of the observed Ha and HI RCs to perform 
the subsequent analysis. 
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Figure 8. Left: NFW fit of Ha DM RC of UGC 4284. Right: NFW fit of Ha DM RC of UGC 11861. In both cases the contours 
contain less than four points. For this reason we do not show that information for these solutions. The number of degrees of freedom 
(DOF) are 78 for UGC 4284 and 43 for UGC 11861. 


halo accomplish better than the DC DM halo. Considering 
the range of disk stellar masses analysed in this study (see 
section and DM halo masses obtained throughout this 
analysis for the DC DM halo our results reproduced closely 
the results presented by Di Cintio et al. (20141 (see figure 4 
of the quoted article), in fact the inner slope outcome from 


the DM RC fitted in this work with the DC DM halo are 
inside the interval [—0.6, —1.2], therefore the results of the 
fits performed with the DC DM halos give in all cases a 
cuspy DM profile. This fact could in part justify the higher 
DM halo masses, with respect to BK DM halo, needed to 
reproduce the DM RCs considered in this work using the 
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R{kpc) R{kpc) 


Figure 9. Left-. BKH fit of Ho DM RC of UGC 4284. Right: BKH fit of Ho DM RC of UGC 11861. The number of degrees of freedom 
(DOF) are 78 for UGC 4284 and 43 for UGC 11861. The contours are displayed in figure [l^ 
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Figure 10. Left: DCH fit of Ho DM RC of UGC 4284. Right: DCH fit of Ho DM RC of UGC 11861. The number of degrees of freedom 
(DOF) are 78 for UGC 4284 and 43 for UGC 11861, the corresponding contours are displayed in figure [T^ 


DC DM halo. In general the fitting errors are about 0.01% 
for the DM halos radii and 0.1% for the DM halos masses, 
in the case where the DM halos fits succeeded. 


In Appendix A we additionally test the Einasto 
|Einasto|l965 1 DM halo (HE) and the Stadel (Stadel et al. 


2009) DM halo (HSD) as feasible options to model the Ha 


and HI DM RCs of UGC 4284 and UGC 11861. Our purpose 
is to compare the outcome from these two DM halos with 
the other results encountered in this research to establish if 
the analysed RCs are best reproduced by cuspy or core pro¬ 
files. As mentioned in Appendix A we regarded the HE and 


HSD DM halos as two free parameters DM halos fixing the 
third parameters for both halos in order to have the same 
degrees of freedom of the other DM halos considered in this 
work. The results showed in the tables Al and A2 clearly 
determine that the cored HE and HSD DM halos fits better 
than BK DM halo and also than DC DM halo practically in 
all cases, with the exception of the DM Ha RC of UGC 4284 
where the DC DM halo produces better results. In summary 
our main result is that the cored DM profiles HE and HSD 
better represent the DM content of UGC 4284 and UGC 
11861. 
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Figure 11. Left: NFW fit of HI DM RC of UGC 4284. Right: NFW fit of HI DM RC of UGC 11861. In both cases the contours 
contain less than four points for this reason we do not show that information for these solutions. The number of degrees of freedom 
(DOF) are 22 for UGC 4284 and 10 for UGC 11861. 
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Figure 12. Left: BKH fit of HI DM RC of UGC 4284. Right: BKH fit of HI DM RC of UGC 11861. The number of degrees of freedom 
(DOF) are 22 for UGC 4284 and 10 for UGC 11861. The contours are displayed in figure [T^ 


7 DISCUSSION 

The ACDM (A Cold Dark Matter) paradigm, predicts cuspy 
density DM profiles in the centre of galaxies, whereas obser¬ 
vational tests, using RCs analysis of spirals, irregular, dwarf, 
and low surface brightness galaxies found that cored density 
DM profiles better represent the central part of observed 
galaxies. The disagreement between the ACDM scenario and 
the observational results, led to consider possible modifi¬ 
cations to the ACDM scheme. One possible refinement of 
the ACDM paradigm consists on the incorporation of bary- 


onic physics on small scales (< 1 kpc). The basic idea to 
transform a cuspy central density profile to a cored one, is 
to remove baryonic matter from the central kpc of galaxies 
through some perturbation to their gravitational potential. 
Supernova explosions and driven-gas outflows, for instance, 
could be very efficient processes to obtain a central DM core 
in the centre of galaxies. Other possible mechanisms are gas 
and metals cooling and star formation. Some numerical and 
analytical attempts have been successful to achieve the de- 
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Figure 13. Left: DCH fit of HI DM RC of UGC 4284. Right: DCH fit of HI DM RC of UGC 11861. The number of degrees of freedom 
(DOF) are 22 for UGC 4284 and 10 for UGC 11861. The contours are displayed in figure 


Chi Square Filled Contours Ha UGC 4284 Burkert DM Halo 



Chi Square Filled Contours HI UGC 4284 Burkert DM Halo 



Figure 14. filled contours at 70% confidence level. First Panel: Ha DM RC of UGC 4284 BKH fit. Second Panel: HI DM RC of UGC 
4284 BKH fit. Third Panel: Ha DM RC of UGC 11861 BKH fit. Fourth Panel: HI DM RC of UGC 11861 BKH fit. The blue crosses 
represent the solutions found through the fitting procedure. 
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Chi Square Filled Contours Ha UGC 4284 DiCintio DM Halo 
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Figure 15. filled contours at 70% confidence level. First Panel. Ha DM RC of UGC 4284 DCH fit. Second Panel HI DM RC of UGC 
4284 DCH fit. Third Panel Ha DM RC of UGC 11861 DCH fit. Fourth Panel HI DM RC of UGC 11861 DCH fit. The blue crosses 
represent the solutions found through the fitting procedure. 


Table 3. Ha DM RCs fit result for UGC 4284 (DOF=78). 


UGC 4284 NFW 

Md + 0'^ 

Md - 30% 

Md + 30%. 

Mh(MQ) 

7.07 X 10^® 

3.75 X IQi® 

3.98 X lOi® 

RjUkpc) 

97.90 

76.33 

65.26 

^Ur-i 

15.36 

20.74 

13.46 

Edfij^ 

1.45 X 10® 

1.30 X 10® 

1.70 X 10® 

UGC 4284 BKH 

Md + 0% 

Md - 30% 

Md + 30%. 

Mh{MQ) 

3.63 X 10® 

2.91 X 10® 

4.52 X 10® 

Rh (kpc) 

10.25 

10.13 

10.45 

Xr 

8.94 

12.72 

5.99 

Edm 

1.22 X 10-1® 

3.04 X 10-1® 

1.74 X 10-1® 

UGC 4284 DCH 

Md + 0% 

Md - 30% 

Md + 30%. 

Mh{MQ) 

5.37 X loll 

7.86 X IQii 

3.57 X loll 

Rh{kpc) 

29.12 

36.15 

22.53 

Xr 

3.32 

2.94 

3.66 

Edm 

3.83 X 10-1® 

7.09 X 10-1-1 

1.48 X 10-1® 

Md = 7.92 X 10®Mq 
^ Reduced x® 

^ Estimated distance from the minimum 

sired DM cored central profile (Governato et al.|2010112012| 


Regardless of all the considerable observational and 


theoretical endeavours to contribute to the solution of the 
cuspy/core controversy, the debate is still open and in this 
work we attempt to discriminate between cuspy and core 
DM halos through DM RCs fitting after subtracting from 
the observed RCs the stellar and gas rotational components. 
We consider high resolution Ha and HI RCs of two irregu¬ 
lar galaxies, UGC 4284 and UGC 11861, available in the 
literature, to examine the mass content of both objects by 
subtracting to the observed Ha the stellar rotation velocity 
and to the HI RCs both the stellar and gas rotation velocity. 
We obtained the stellar rotational component through stel¬ 
lar population synthesis of the disk of these two galaxies. 
We recognise that in the case of both galaxies, once sub¬ 
tracted the determined stellar and gas rotation velocity, it 
is missing more than 50% of mass, for a stellar disk mass of 
Md = 7.92 X 10® Mq in the case of UGC 4284 and more than 
40% of mass for a stellar disk mass of Md = 5.38 x 1O^®M0 
in the case of UGC 11861. We explored the possibility that 
DM could account for the missing mass in these galaxies, 
and we fitted the derived DM RCs with five DM models, 
namely the NFW DM halo, the BK DM halo, the DC DM 
halo, the HE DM halo and the HSD DM halo, to estimate 
the assumed DM contribution to the total mass of UGC 
4284 and UGC 11861. 

The cuspy DC DM halo fit better the Ha DM RC of 
UGC 4284 than any other DM halos used in this study, in 
all the other cases, say HI RC of UGC 4284 and Ha and HI 
RCs of UGC 11861, the cored DM models give a better fit to 
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Table 4. Ha DM RCs fit result for UGC 11861 (DOF=43). 


Table 6. HI DM RCs fit result for UGC 11861 (DOF=10). 


UGC 11861 NEW 

Md + 09i 

a 

Md - 30% 

Md + 30%. 

MhiMQ) 

2.93 X 10 

12 

7.05 X 10^® 

4.79 X lO^® 

RfMpc) 

54.26 


97.91 

65.27 


70.86 


75.67 

55.67 

Edl^ 

9.62 X 10* 

1.37 X 10® 

5.58 X lO^® 

UGC 11861 BKH 

Md + 0% 

Md - 30% 

Md + 30%. 

Mh{MQ) 

3.16 X 10® 

3.25 X 10® 

3.08 X 10® 

Rh (kpc) 

9.16 


10.21 

8.29 

Xr 

46.10 


58.32 

35.71 

Edm 

1.40 X 10“ 

3.01 X 10“^® 

9.10 X 10-1® 

UGC 11861 DCH 

Md + 0% 

Md - 30% 

Md + 30%. 

Mh(MQ) 

2.89 X 10 

12 

4.09 X 10^® 

1.96 X lOi® 

Rh{kpc) 

71.24 


88.29 

54.92 

Xr 

32.08 


42.26 

25.65 

Edm 

2.40 X 10“ 

1.40 X 10-1"^ 

6.03 X 10-1® 


^ Md = 5.38 X 
^ Reduced 

^ Estimated distance from the minimum 


Table 5. HI DM RCs fit result for UGC 4284 (DOF=22). 


UGC 4284 NEW 

Md + 0% 

a 

Md - 30% 

Md + 30%. 

Mh{MQ) 

4.16 X 10 

12 

3.64 X lOi® 

2.94 X IQi® 

RiJkpc) 

97.90 


97.90 

79.03 


305.74 


299.16 

307.61 

eJI^ 

4.77 X 10® 

4.52 X 10* 

2.61 X 10® 

UGC 4284 BKH 

Md + 0% 

Md - 30% 

Md + 30%. 

Mh{MQ) 

1.0 X 10 

s 

1.0 X 10* 

1.0 X 10* 

Rh (kpc) 

1.73 


1.85 

1.62 

Xr 

108.25 


131.65 

118.67 

Edm 

1.34 X 10 

2.68 X lO-l” 

1.67 X 10-* 

UGC 4284 DCH 

Md + 0% 

Md - 30% 

Md + 30%. 

Mh{MQ) 

1.13 X 10 

11 

1.26 X lOii 

9.83 X IQi® 

Rh(kpc) 

25.01 


31.69 

19.08 

Xr 

82.96 


103.41 

77.54 

Edm 

3.18 X 10" 

-11 5.67x10-11 

2.67 X 10-11 


UGC 11861 NEW 

Md + 07i 

a 

Md - 30% 

Md + 30%. 

MhiMe) 

5.46 X 10 

12 

6.24 X IQi® 

3.89 X lOi® 

Rhjkpc) 

81.36 


88.45 

78.67 


97.35 


95.34 

92.67 

Ed50 

3.28 X 10* 

5.32 X 10* 

6.54 X 10® 

UGC 11861 BKH 

Md + 0% 

Md - 30% 

Md + 30%. 

Mh{MQ) 

6.53 X 10* 

6.12 X 10* 

6.94 X 10® 

Rh (kpc) 

3.87 


4.01 

3.75 

Xr 

37.60 


38.35 

36.93 

Edm 

2.13 X 10” 

■1® 2.53 X 10-1® 

1.89 X 10-1® 

UGC 11861 DCH 

Md + 0% 

Md - 30% 

Md + 30%. 

MhiMe) 

1.35 X 10 

12 

1.70 X IQi® 

1.01 X lOi® 

Rh{kpc) 

68.38 


85.26 

52.22 

Xr 

44.98 


54.95 

48.08 

Edm 

4.17 X 10” 

-1® 3.59 X 10-1® 

6.15 X 10-1® 


^ Md = 5.38 X lO^^M© 

^ Reduced x^ 

^ Estimated distance from the minimum 


As mentioned in section [^the cuspy result we obtain in the 
case of the Ha DM RC of UGC 4284 is in accordance with 
the numerical study of |Di Cintio et "M? (20141 (see figure 
4 of that article), given that our range of stellar, DM halo 
masses and inner slopes coincide with those presented in the 
quoted hgure. Another important outcome of this study is 
the fact that the M/L ratio is not constant in none of the 
two studied galaxies (see figures and |^ . 

At this point we can answer, at least for UGC 4284 
and UGC 11861, to the question we asked in the Introduc¬ 
tion, i.e., if the DM prohle slope could change substantially 
obtaining the stellar disk from stellar population studies in¬ 
stead of assuming a constant M/L as previous studies done 
in the past (see Introduction and references therein). From 
the results outlined above it seems that in the case of UGC 
4284 and UGC 11861, the answer to this question is a nega¬ 
tive one. In the future we plan to repeat this study for other 
high resolution RCs, in order to asses this question for a 
larger number of galaxies. 


Md = 7.92 X 10®Mq 
^ Reduced y® 

Estimated distance from the minimum 

the DM RCs of both galaxies than the cuspy DM profiles. 
In more detail the HE and HSD DM halos with two free 
parameters perform better than the BK DM halos except 
for the Ha DM RC of UGC 4284 with stellar disk mass 
Md + 30%. The HE DM halo achieve better results than the 
HSD DM halo for the HI DM RCs of both galaxies and also 
for the DM Ha RC of UGC 11861 with stellar disk mass 
Md + 0% and Md -I- 30%. The HSD DM halo give better 
results than the HE for the DM Ha RC of UGC 4284 and 
DM RC of UGC 11861 with stellar disk mass Md — 30%. 


8 CONCLUSIONS 

In this study we analyse the inner and outer mass distribu¬ 
tion of the two irregular galaxies UGC 4284 and UGC 11861, 
by means of high resolution Ha and HI RCs available in the 
literature. Additionally, we constrained the stellar content 
of UGC 4284 and UGC 11861 throughout stellar population 
synthesis studies, applied to the photometric observations of 
these two galaxies in several Hlters, we also derive HI+He 
RCs for the HI observed RCs of both galaxies and we sub¬ 
tracted the derived stellar disk from the observed Ha RCs 
and the stellar and gas disk from the observed HI RCs. 

We fitted the determined DM RCs for both galaxies, 
examining the NEW, BK, DC, HE and HSD DM models. 
The principal results of this research can be summarized 
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according to the following points: 

(i) Cored DM models better fit the observed shape of 
the two DM Hq and HI RCs considered in this study, with 
respect to cuspy DM models, except for the Ha DM RC of 
UGC 4284 where the cuspy DC DM halo is the preferable 
model. 

(ii) Cored exponential DM model Einasto reproduced 
better the DM HI RCs of UGC 4284, and UGC 11861, than 
the other models tested in this paper. 

(iii) Cored exponential DM model Einasto reproduced 
better than all the other models the DM Ha RCs of UGC 
11861 for stellar masses Md + 0% and Md + 30%. 

(iv) The cored Stadel DM model, in the case of the Ha 
DM RC of UGC 4284, performed better than Einasto DM 
models for all the stellar masses considered. 

(v) The cored Stadel DM model, in the case of the Ha 
DM RC of UGC 11861 and a stellar mass Md — 30%, 
achieved better results than Einasto DM models. 
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APPENDIX A: EINASTO AND STADEL DM 
RC FIT RESULTS 


We also consider the Einasto (HE) | Einasto]|1965 1 and the 
Stadel ( Stadel et al.|2009 1 (HSD) DM velocity profiles in our 
endeavour to reproduce the Ha and HI DM RCs of UGG 
4284 and UGG 11861. The motivation to employ the HE 
and HSD DM halos are emphasized below, even if at the 
moment it does not exists a physical justification of these 
two DM profiles. During the fitting process we fixed the 
third parameter of both DM halos, therefore the number 
of degrees of freedom is the same of the other DM halos 
examined in this work. 

The HE was introduced for the first time by |Einasto| 
(19651 to build a model of the distribution of stars in the 
Galaxy. Navarro et al. (2004 20101 established that the HE 
fits the DM density profiles produced in the framework of 
AGDM cosmological numerical simulations better than the 
NFW profile. Ghemin et al. (20111 used for the first time the 
HE DM models to fit 34 nearby galaxies, determining that 
the HE DM halo represented better the observed RCs of the 
studied galaxies than the Pseudo-Isothermal DM model. 

The HE represents a three parameters model in the form 
of an exponential-like profile and also incorporates cusped 
and cored DM halos. We used the approximation of ICiottil 


& Bertin (19991 for the dimensionless scale factor ci„ 

The HSD have been proposed by [Stadel et al.| ( |2009| ) 
accomplishing several AGDM Galactic mass DM halo nu¬ 
merical simulations, the authors concluded that the HSD 
better reproduced the shape of the simulated DM halos den¬ 
sity profiles also if compared with the HE. The HSD is an 
exponential-like profile with three free parameters, but ac¬ 
cording to [Stadel et al.j ( |2009[ ) can also be regarded as a 
DM density profile with two free parameters with functional 
form: 


p(r) 




where R\ is the HSD scale radius and A is a parameter 
that governs the shape of the HSD. In this study we also 
provide an exact derivation of the HSD velocity profile in 
Appendix B. In particular the relation (BIO) of Appendix 


B, is a tabulated integral from Gradshteyn & Ryzhik 120071. 


The relation (B19) of the same Appendix B, is the expression 
we used in this investigation to fit the DM RCs of UGG 4284 
and UGG 11861 employing the HSD. The results of the HE 
and HSD fits are listed in table Al and table A2. 


APPENDIX B: CIRCULAR VELOCITY 
DERIVATION FOR THE STADEL DM MODEL 


The spatial density relation for the HSD halo density profile 
is expressed by the relation: 

p{r) = Poe . (Bl) 

We can write this spatial density in the following form: 


pR) = Po 




(B2) 


© 2002 RAS, MNRAS 000 , [^?? 









































Outer and inner mass distributions of UGC 4284 and UGC 11861 


17 


Table Al. Ho DM RC fit results (UGC 4284 DOF=78 AND UGC 11861 DOF=43). 


UGC 4284 EINASTO 

Mi + 0% 

Mi - 30% 

Mi + 30%. 

MhiMe) 

1.22 X 10^2 

1.26 X 10i2 

1.22 X 10i3 

Rh{kpc) 

100 

100 

100 

rih 

3.16 

2.92 

3.37 

Xr 

4.43 

4.04 

4.96 

Edm 

8.37 X 10“® 

2.40 X lO-l” 

4.10 X lO-l" 

UGC 4284 STADEL 

Mi + 0% 

Mi - 30% 

Mi + 30%. 

Mh{MQ) 

8.70 X lO^i 

1.14 X 10i3 

5.17 X lOii 

Rh(kpc) 

47.13 

59.65 

33.73 


3.79 

3.80 

3.81 

Xr 

6.13 

7.10 

8.77 

Edm 

1.07 X 10-11 

1.61 X 10-1-1 

6.39 X 10-® 

UGC 11861 EINASTO 

Mi + 0% 

Mi - 30% 

Mi + 30%. 

Mh{MQ) 

1.19 X lOii 

8.27 X IQil 

1.81 X lOii 

Rh(kpc) 

14.92 

12.85 

18.12 

rih 

1.13 

0.92 

1.39 

Xr 

6.18 

10.04 

3.74 

Edm 

3.81 X 10-1"! 

1.68 X 10-13 

3.88 X 10-1® 

UGC 11861 STADEL 

Mi + 0% 

Mi - 30% 

Mi + 30%. 

Mh{MQ) 

5.27 X lOii 

4.61 X IQii 

9.86 X loll 

Rh(kpc) 

31.23 

31.22 

41.89 

Ah 

3.65 

3.59 

3.75 

Xr 

8.05 

8.02 

8.96 

Edm 

1.74 X 10-13 

2.70 X 10-13 

4.11 X 10-3 


Table A2. HI DM RC fit results (UGC 4284 DOF=22 AND UGC 11861 DOF=10). 


UGC 4284 EINASTO 

Mi + 0% 

Mi - 30% 

Mi + 30%. 

Mh(M0) 

1.82 X lOio 

1.59 X IQi® 

2.06 X lOio 

Rh(kpc) 

5.42 

5.23 

5.64 

rih 

0.80 

0.73 

0.87 

Xr 

5.90 

5.74 

6.10 

Edm 

1.68 X 10-13 

1.59 X 10-13 

5.99 X 10-13 

UGC 4284 STADEL 

Mi + 0% 

Mi - 30% 

Mi + 30%. 

MhiMe) 

7.41 X 10i° 

6.44 X IQi® 

8.02 X lOio 

Rh{kpc) 

12.18 

11.60 

11.84 

Ah 

3.82 

3.77 

3.73 

Xr 

25.30 

24.33 

17.95 

Edm 

3.52 X 10-13 

4.78 X 10-13 

1.68 X 10-1"! 

UGC 11861 EINASTO 

Mi + 0% 

Mi - 30% 

Mi + 30%. 

Mh(M0) 

5.70 X lOio 

5.16 X IQi® 

6.54 X lOio 

Rh(kpc) 

9.20 

10.30 

11.40 

rih 

1.12 

1.10 

1.15 

Xr 

19.50 

17.25 

20.51 

Edm 

1.26 X 10-3 

1.30 X 10-3 

1.17 X 10-3 

UGC 11861 STADEL 

Mi + 0% 

Mi - 30% 

Mi + 30%. 

Mh(M0) 

3.65 X loll 

3.31 X IQii 

4.0 X IQii 

Rh{kpc) 

26.60 

27.04 

26.18 

Ah 

3.79 

3.79 

3.79 

Xr 

26.70 

31.38 

27.98 

Edm 

2.06 X 10-1"! 

1.44 X 10-13 

1.21 X 10-1^ 
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using exponential properties we obtain the relation: 

1 

1 

Pp) = po 




(B3) 


thus we can use this form of the spatial density to obtain 
the mass density up to a radius r, with the relation: 


rr 

M{r) -A-k p{s)s^ds 
Jo 

substituting in this relation we obtain: 


(B4) 


M{r) = 47rpo 


rr 

Jo 






ds (B5) 


^ = 


in s = R\ (s = i?A ^6”^ “ 8'iid 


with the substitution e“ = 
we obtain s = 

u = -Aln(^l + => du = -X-fjn (l + ^ du = 

1 

—A , —V ds ^ du = —-PeT^ ds ^ ds = du. 

(l+^) «A A 


For the limits of integration we have that u = In 


we obtain w = 0 for r = 0 and m = In 




(^+rfy) 

for r = r. 


It is important to note that the upper integration limit 


(' + *)■ 


is negative given that u = In 

—Ain with {\,r ,Rx) > 0. 

We proceed in the calculation of the integral substituting all 
the derived quantities in the above relation obtaining: 


M(r) = - 


47rpoi?i 


A 


MU<V) ^2 y g 


e ^ ds (B6) 


After simplifying this relation and making the square we 
obtain the following relation: 


M^r) = r e-^ (e-^P + - 26“V) ds 

^ J{u<0) ^ ' 

(B7) 

multiplying and simplifying we get the expression: 


M{r) = 


47rpo-Ri 


/ 

d (l 


(u<0) 


^+3£ 

A 


+ e 


- 2e 


“^+25 

A 


(B8) 

naming the first integral in the above expression as A(s) = 


r 

J(u<0) 


ds, the second B{s) = /(° <o) s 


ds 


and the third integral C{s) = e a J ds we can 

rewrite our mass integral in the following manner: 


At this point considering each integral separately and con¬ 
sidering the indefinite integral given by the expression: 


f e-dx = exp (] erf y/dx + ^ 

J 2\ a \ a J I xfd 

(BIO) 

where the erf{x) is given by relation erf(x) = 

in the case of the hrst integral A(s), we obtain a = ^; 
& = For the second integral B(s) in a similar way we get 
a = ^; fe = ^ and for the third integral we obtain a = ^; 
b = J. All the integrals have c = 0. 

Thus the previously indefinite integrals A(s), B{s) and 
C{s) transform according to the relation: 


A(s) = I < 




B{s) = Je V ^ )ds = exp ) erf 


C{s) = y e ( ^ )ds= exp ) erf 


1 3 

A^^27a 

(Bll) 


A''+2yA 

(B12) 


1 1 

H- 'p- 

A \/A 

(B13)' 


now, we calculate the definite integrals setting u = —v 
where u > 0 with the reminder that erf{—x) = —erf{x) we 
obtain the following results: 


f 

J(-v<0) 


e ^ ^ ) ds = A(0) — A{—v) = exp ^ 


4A 


r 

fl 3 


r 3 ] 

erf 

\ —V -= 

+ erf 

.2VA. 

V A 2VX 


/: 


0 

(—■i;<0) 


e A ^ ) ds = B(0) — B{—v) = exp 


X < erf 


a"" 2^ 


-I- erf 


1 

“4A 

1 

2^. 


f 

J(-v<0) 


e V ^ ) ds = (7(0) — C(—v) = exp f—^ 

2 \ A 


ds 


X < erf 




-I- erf 


VaJ 


Let us define the part involving the erf functions in the 
above expressions in the following way: 


hrfiiv, A) = erf 


a’' 2yA 


-I- erf 


2yA. 


(B14) 


M(r) = 


dTrpod^A 

A 


[A(s) + B{s) - 2(7(s)] (B9) 


hrfoiv, A) = < erf 


a"" 2Va 


-I- erf 


2yA. 


(B15) 
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lerfsG, 0 = I erf 

at this point we can rewrite the HSD mass function in 
the following form: 

M(r) = [e-i^hrfi{v,\) + e-^hrf2{v,X)-2e-^hrf3iv,\)\ 

(B17) 

where —v = In ^^ and the circular velocity is 
given by V'c(r) = GM(r)/r so the above relation becomes: 

[e-A7e./i(w,A) + e-K7^,^2(v,A)-2e-i7e./3(w,A)] 

(B18) 

we can express the cicular velocity expression as a func¬ 
tion of the initial mass of the system remembering that 
M\ = poVa, where Va = \t^R\ is the volume of a sphere 
of radius r containing the mass M for definition, so we ob¬ 
tain: 

Vf{r) = [e-i^7e./i(u.A) + e-^7e./2(u,A)-2e-i7e./3(^;,A)] 

(B19) 

where Mx = ^ivR^po is the mass enclosing the initial 
density po inside a radius R\ and G is the gravitational con¬ 
stant expressed in astrophysical units kpc krrds~^). 




-I- erf 


Va 


(B16) 
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